Mean-field dynamics of a Bose-Einstein condensate in a time-dependent triple-well 
trap: Nonlinear eigenstates, Landau-Zener models and STIRAP 



in 
o 
o 

> 

o 
m 

> 
in 
oo 

t> 
o 
in 
o 

^ : 

Qk 
^— > ■ 
G ■ 
ctf 

3 : 



E. M. Graefe, H. J. Korsch, and D. Witthauf] 

Technische Universitat Kaiserslautern, FB Physik, D-67653 Kaiserslautern, Germany 

(Dated: February 1, 2008) 

We investigate the dynamics of a Bose-Einstein condensate (BEC) in a triple-well trap in a 
three-level approximation. The inter-atomic interactions are taken into account in a mean-field 
approximation (Gross-Pitaevskii equation) , leading to a nonlinear three-level model. New eigenstates 
emerge due to the nonlinearity, depending on the system parameters. Adiabaticity breaks down if 
such a nonlinear eigenstate disappears when the parameters are varied. The dynamical implications 
of this loss of adiabaticity are analyzed for two important special cases: A three level Landau- 
Zener model and the STIRAP scheme. We discuss the emergence of looped levels for an equal- 
slope Landau-Zener model. The Zener tunneling probability does not tend to zero in the adiabatic 
limit and shows pronounced oscillations as a function of the velocity of the parameter variation. 
Furthermore we generalize the STIRAP scheme for adiabatic coherent population transfer between 
atomic states to the nonlinear case. It is shown that STIRAP breaks down if the nonlinearity exceeds 
the detuning. 
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INTRODUCTION 



The experimental progress in controlling Bose-Einstein 
condensates (BECs) has led to a variety of spectacular re- 
sults in the last few years. At very low temperatures, 
the dynamics of a BEC can be described in a mean- 
field approximation by the Gross-Pitaevskii (GPE) or 
nonlinear Schrodinger equation (NLSE) pj. Previously, 
several authors investigated the dynamics of the NLSE 
for a double- well potential in a two-mode approximation 
HSIlllUHSSI- Novel features were found e.g. the 
emergence of new nonlinear stationarystates Q and a 
variety of new crossing scenarios (cf. Q and references 
therein) . Studies of the quantum dynamics beyond mean- 
field theory were reported in 0, OH • First approaches 
of the coherent control of BECs in driven double-well 
potentials have been reported [|J. Another relevant ap- 
plication of the two-level NLSE is the dynamics of a BEC 
in an accelerated or tilted optical lattice PJ || . Further- 
more, the NLSE describes the propagation of light pulses 
in nonlinear media and its two-level analogon has 
been studied in the context of a nonlinear optical direc- 
tional coupler This equation is also known as the 
discrete self-trapping equation and has been applied, e.g., 
to the dynamics of quantum dimers. In fact, the char- 
acteristic loop structures [5( important in the following, 
have been already observed in this context flU ]. 

In the present paper we will extend these studies to 
nonlinear thrcc-lcvcl quantum systems with special re- 
spect to the breakdown of adiabatic evolution due to the 
nonlinearity. In fact, we investigate the dynamics of a 



BEC in a three-level system 



(1) 



in a mean-field approach. The dynamics of the coeffi- 
cients is described by the discrete NLSE 



H(\a\\U 




(2) 



with the nonlinear three-level Hamiltonian 
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H(\af,\b\ 2 ,\c\ 2 ) = 



e + g\a\ v 
v g\b\ 2 







w | . (3) 

to <5 + g|c| 2 
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The dynamics preserves the normalization, which is fixed 
as \a\ 2 + \b\ 2 + |c| 2 = 1. Throughout this paper units are 
rescaled such that all properties become dimensionless 
and h = 1. 

An experimental setup leading very naturally to the 
Hamiltonian is the dynamics of a BEC in a triple- 
well potential. In this case the basis states ip n (x) are 
localized in the three wells. Then e and S are the on- 
site energies of the outer wells and v and w denote the 
tunneling matrix elements between the wells. These pa- 
rameters can be varied by controlling the depth or the 
separation of the wells. The energy scale is chosen such 
that the on-site energy of the middle well is zero. This 
situation is illustrated in figure ^ 

A detailed discussion of these approximations and the 
validity of the model can be found in Q for the case of 
a double-well potential. Furthermore, the existence of a 
dark state for ultracold atoms and molecules, which was 
demonstrated experimentally very recently, could be de- 
scribed successfully within a nonlinear three-level model 
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FIG. 1: Illustration of the triple-well potential and the cor- 
responding three-level model (JSJl investigated in the present 
paper. 



(see [l4j and references therein) Further applications in- 
clude the dynamics of three-mode systems of nonlinear 
optics or the analysis of the excitonic-vibronic coupled 
quantum trimer. 

As previously shown for the two-level system, the non- 
linearity leads to the emergence of new eigenstates with- 
out a linear counterpart, loop structures and novel level 
crossing scenarios. The conce pt o f adiabaticity is very 
different in the nonlinear case |15|. leading to nonlinear 
Zener tunneling JS\ and possibly to dynamical instability 
||. These issues are discussed in section lTTll for the three- 
level system @. Furthermore we analyze adiabatic co- 
herent population transfer. In the linear case a complete 
population transfer can be achieved using the Stimulated 
Raman Adiabatic Passage (STIRAP) via a dark state of 
the system ^||. Coherent control techniques for ultra- 
cold atoms in a triple-well trap were previously investi- 
gated by Eckert et. al. for the case of single atoms ]vf\ - 
In the present paper it is shown that a mean-field inter- 
action plays a crucial role and that the STIRAP scheme 
fails if the nonlinearity exceeds a critical value due to the 
breakdown of adiabaticity. 



II. NONLINEAR EIGENSTATES 

Nonlinear eigenstates and eigenvalues can be defined 
as the solutions of the time-independent NLSE 



H(\a 




(4) 



with the chemical potential fi. Obviously an interpreta- 
tion in terms of linear algebra is not feasible anymore. 
Therefore - although equation (0J looks like a common 
eigenvalue-equation - most of the implications of linear 
quantum mechanics are not valid anymore. Nevertheless 
these eigenstates are of great importance for understand- 
ing the dynamics of the system, as they still are station- 
ary solutions of the time dependent NLSE J2J • 



In the two-level case one can calculate the eigenstates 
by solving a forth-order polynomial equation 5]. The 
case of three levels turns out to be a bit more compli- 
cated. As all parameters are real, the amplitudes (a, 6, c) 
are also real. For almost all parameters (except in the 
"uncoupling" limits v, w — > or 5, e — > oo) the ampli- 
tudes are all non-zero. Thus the variables x = b/a and 
y = c/b are well defined and determined by the equations 

(1 - x 2 y 2 ) (g + 5 + e + v(x + x^ 1 ) + w(y + y" 1 )) 
-3x 2 +y 2 (e + vx)-3wy- 1 + 35 = (5) 

(1 - x 2 ) (g + S + e + v(x + x^ 1 ) + w(y + y^ 1 )) 
-3a;" 1 -3wy + 3x 2 (e + vx) = 0. (6) 

Equation © can be solved for y explicitly. Substitution 
of the result into equation (JSJ then leads to a single equa- 
tion for x, which can be solved numerically. In the limit 
of large S and e, the eigenvalues are easily found by non- 
linear optimization with the linear eigenvalues as initial 
guesses. 

Similar to the case of two levels Q, the nonlinear 
three-level model can also be described as a classical 
Hamiltonian system. Introducing the variables p\ = \a\ 2 , 
P3 = |c| 2 , qi = arg(6) - arg(a) and q 3 = arg(6) - arg(c), 
the dynamics is given by the conjugate equations 



Pi 



on 

= — — — and 



1j = 



&H 

dp 3 



(7) 



with the classical Hamiltonian function 

H = epi + S P3 + | (p? + p 2 3 + (1 - Pl - P3 ) 2 ) (8) 



+2^1 -pi-p 3 {v^/p : [cos(q 1 ) + w % /pi'cos(g 3 )) 

and the normalization condition \b\ 2 = 1 — \a\ 2 — \c\ 2 . 
The eigenstates (@J, resp. the stationary states of the 
system, are given by <jj = P j = 0. Hence they correspond 
to the critical points of this classical Hamiltonian, which 
one finds to be given by the real roots of two polynomials 
in pi and p 3l one of them of 8th order in p\ and of 7th 
in p 3 the other one the opposite way round. At these 
critical points, the classical Hamiltonian and the chemical 
potential jx are simply related by 



n = v-^(pi+pj + (i-pi- P 3) 2 ) 



(9) 



New eigenstates emerge when the nonlinearity \g\ is 
increased, starting from three eigenstates in the linear 
case g = 0. From the point of view of standard quan- 
tum mechanics, the existence of additional eigenstates is 
a strange issue. Simply considering them as stationary 
states is a more adequate point of view. One can in- 
tuitively understand that a strong interaction compared 
to the on-site-energies not only modifies the stationary 
states of the linear case but can stabilize other states as 
well and thus create additional eigenstates. For a deeper 
insight in the features of these states see, e.g., 0, where 



3 



the nonlinear eigenstates in a symmetric double- well trap 
are studied in detail. 

In terms of the corresponding Hamiltonian system no 
conflict arises at all. The (dis)appearance of critical 
points due to a variation of the system parameters is 
a well-known issue in the theory of classical dynamical 
systems. In fact this is a major ingredient of catastrophe 
theory [Isj . Another important results from the theory 
of classical dynamical systems is that the number of ex- 
trema minus the number of saddle points of Ti is constant 
(see, e.g. Thus two fixed points, one elliptic and one 
hyperbolic, always emerge together. 

Furthermore this classical description is especially 
suited for the discussion of adiabatic processes [15j ■ Note 
that the dynamics of the discrete nonlinear Schrodingcr 
equation can become classically chaotic [2Cj . 

III. LANDAU-ZENER TUNNELING 

In this section we study the evolution of the system 
under variation of the parameters in the manner of the 
Landau-Zener model, considering v and w independent of 
time. In particular we will focus on the so called equal- 
slope case [2l| , for which S is constant in time as well and 
e = at. In the linear case g = the instantaneous eigen- 
values in dependence of e = at (the so called adiabatic 
levels) form two consecutive avoided crossings which ap- 
pear at e = A± with gaps of size = —^i,2v/ri2,i with 

Ai, 2 = 8/2±yj8 2 /A + w 2 and m, 2 = W / A| 1 + w 2 . An ex- 
ample is shown in figure [21 (dashed lines). If the system 
is prepared in a state on the lowest branch for t — ► — oo 
and the parameters are varied infinitely slow (i.e. a — > 0) 
with g — 0, the adiabatic theorem states that the quan- 
tum state will follow the adiabatic eigenstate up to a 
global phase . However, for a finite value of a, parts 
of the population will tunnel to the other adiabatic lev- 
els, mainly at the avoided crossings. In the case of the 
linear two-level system the transition probability P be- 
tween the two adiabatic levels is given by the celebrated 
Landau-Zener formula [2^| 

Plz = exp (~2nv 2 /a) . (10) 

For the equal-slope case in linear three-level systems, one 
can show that the transition probability to the highest 
adiabatic eigenstate (or equivalently the survival prob- 
ability in the diabatic eigenstate populated initially) is 
also given by the Landau-Zener-formul a 1J10I) and con- 
sequently independent of w and S (see [21], |2J, |25j and 
references therein). 

In the following we study the behavior of the system 
in the equivalent scenario with g ^ 0. For that purpose 
we first calculate the nonlinear eigenvalues, which are de- 
fined by equation As in the two- level case loops 
emerge with increasing nonlinearity near the points of 
the avoided crossings at critical values of \g\ = g c \ resp. 
\g\= g C 2, which are governed by the width of the gap at 
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FIG. 2: Nonlinear eigenvalues as defined by equation in 
dependence of e for S = —0.4, v = 0.1, w = 0.2 and g = —0.4 
(solid line) resp. g = (dashed line). 
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e = a t 

FIG. 3: Time evolution of the population in the first level 
( | ot(£) | 2 , open circles o— ) in comparison with the first compo- 
nent of the corresponding adiabatic eigenstates (|a a( j| 2 , solid 
line) for e = at with a = 0.001 and 6 = -0.4, v = 0.1, to = 0.2 
and g = —0.4. 

the avoided crossings in the linear case. An estimate of 
the critical values g c \ and g c i is obtained by considering 
the two avoided crossings as independent and hence ne- 
glecting the influence of the remaining third level. For a 
two-level system, the critical nonlinearity can be calcu- 
lated exactly which yields g c \ — 2v\ and g c i = 2v%. 
Our numerical studies show that this is an acceptable ap- 
proximation and at least a lower bound to the real value. 
An example of the nonlinear eigenvalues in dependence of 
e is shown in figure[21for S = —0.4, v = 0.1 and w = 0.2. 
We choose g = —0.4 < such that the loops emerge on 
top of the two lower adiabatic levels for better compar- 
ison with However, similar results are found if the 
signs of both g and a are altered. 

The first component of the nonlinear eigenstates as- 
sociated to the lowest level in figure [3 are illustrated in 
figure^ (solid line) . The loop structure of the eigenvalues 
manifests in an S-shaped structure of the components of 
the corresponding eigenstate ^4] . For a better comparison 
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FIG. 4: Tunneling probability P(a) in dependence of the 
velocity of the parameter variation for g = —0.03 (open cir- 
cles), g = —0.16 (full circles) and g = —0.4 (open squares) 
and 8 = —0.4, v — 0.1, w — 0.2. The solid lines are drawn to 
guide the eye. The Landau-Zener formula 110H for the linear 
case is shown for comparison (dashed line). 



between the figures and we included the labels A-D. 
The nonlinear eigenstates emerge/ vanish in a bifurcation 
at edges of the loops resp. the S-structure (marked by B 
and C in the figures). 

As in the two-level system, the appearance of the loops 
leads to the breakdown of adiabatic evolution. To illus- 
trate this issue we calculate the dynamics for the same 
parameters as above with a slowly varying e = at with 
a = 0.001 and compare them to the relevant instanta- 
neous eigenstates. The system is initially (t — > — oo) 
prepared in an eigenstate corresponding to the lowest 
level in figure [5] (point A) . The resulting dynamics of 
the population in the first level |a(i)| 2 is shown in figure 
01 (open circles). The solid line shows the population in 
the first well |a a( j| 2 for the corresponding instantaneous 
eigenstates. One clearly observes the breakdown of the 
adiabatic evolution when the adiabatic eigenstate van- 
ishes in a bifurcation at the edge of the S-structure at 
the point B around at w —0.25. 

Due to this breakdown of adiabaticity the transition 
probability P does not vanish for a — > for \g\ > g c . This 
is illustrated in fig. 0] where we plotted P in dependence 
of a for different values of the nonlinear parameter g. 
For weak nonlinearities, e.g. g = —0.03 in the figure, P is 
increased in comparison to the linear case, however it still 
tends to zero in the adiabatic limit, i.e. P — > for a — > 0, 
as no loops have occurred yet. This is no longer true after 
the appearance of the two loops shown in fig. [3 such that 
P(a — > 0) > (cf. fig. . These features are well-known 
from the nonlinear two- level model However, a novel 
feature is the appearance of pronounced oscillations of 
the transition probability P{a) for small values of a due 
to the nonlinear interaction between the different levels. 

As mentioned above, the transition probability for the 
linear case is independent of w and 5. This does not hold 
in the nonlinear case any longer, since the influence of the 
nonlinearity, e.g. the emergence and the structure of the 
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FIG. 5: Tunneling probability P(a) in dependence of the 
velocity of the parameter variation for w — 0.2 (open circles), 
w = 0.4 (full circles) and w = 0.6 (open squares) and S = 
—0.4, v = 0.1 and g — —0.4. The solid lines are drawn to 
guide the eye. The Landau-Zener formula 1101 for the linear 
case is shown for comparison (dashed line). 



loops, obviously depends on v, w and S. This is confirmed 
by our numerical results shown in fig. UJ where P(a) is 
plotted for different values of w. Again pronounced os- 
cillations of P(a) are found for small values of a. These 
oscillations and their dependence on the system's param- 
eters will be studied in detail in a subsequent paper. 



IV. NONLINEAR STIRAP 

The STIRAP method, primarily proposed and real- 
ized for atomic three-level systems [T||, allows a robust 
coherent population transfer between quantum states. In 
the meanwhile STIRAP has been generalized to systems 
with multiple levels and the preparation of coherent su- 
perposition states |26( . In this case the coupling between 
the atomic bare states are realized by slightly detuned 
time-dependent laser fields with Rabi frequencies v(t) 
and w(t). In the rotating wave approximation at the 
two-photon resonance, the dynamics of the three-level 
atom is given by the Hamiltonian © with g = and a 
fixed detuning A = — 5 = — e. Using the STIRAP scheme 
one can achieve a complete population transfer from level 
to level ips by an adiabatic passage via a dark state 
^ds of the system, which is a superposition of and ipz 
alone. If the coupling w(t) between the levels tp2 and 
ipz is turned on before the coupling v(t) between the lev- 
els ipi and ip2, the system's dark state is rotated from 
tpds(t = — oo) = ipi to ipds(t = +oo)) = 1P3. If the param- 
eters v(t) and w(t) are varied sufficiently slowly, the sys- 
tem can follow the dark state adiabatically which leads to 
a complete population transfer from level ip\ to ip^. This 
counterintuitive pulse sequence of the coupling elements 
v(t) and w(t) is illustrated in fig. H3 (upper panel). 

However, the situation is more involved in the nonlin- 
ear case. To begin with we consider the case g, A > 0. As 
noted in the previous section, new nonlinear eigenstates 
emerge if the nonlinearity \g\ exceeds a critical value g c 
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FIG. 6: Eigenvalues of the nonlinear Hamiltonian for a 
detuning A = 0.1 and a nonlinearity g — 0.2 (lower panel) for 
time-dependent couplings v(t) and w(t) (upper panel). 
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FIG. 8: Time evolution of the population in the first (|a(t)| 2 , 
open circles) and the third level (|c(t)| 2 , crosses) in compari- 
son with the population in the first level for the instantaneous 
eigenstates (|a a d| 2 , solid lines) for A = 0.1 and g — 0.2 and 
the couplings v(t),w(t) shown in fig. [H] (upper panel). 
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FIG. 7: Left: Magnification of fig. |S| Right: Squared modu- 
lus of the first component (|a a d| 2 ) of the corresponding eigen- 
states. The avoided crossing that leads to a breakdown of 
adiabaticity is marked by an arrow. 



depending on the other sytem parameters. This can give 
rise to new crossing scenarios leading to a breakdown of 
adiabaticity. To illustrate this issue, the eigenvalues and 
eigenstates are calculated for a detuning A = 0.1 and 
g = 0.2. The resulting adiabatic eigenvalues as well as the 
STIRAP pulse sequence of the couplings v(t) and w(t) are 
shown in fig. At a first glance this picture looks quite 
similar to the linear case. The levels are shifted slightly 
due to the mean-field energy and a few additional non- 
linear eigenvalues emerge for large \t\, because then the 
coupling elements v and w are small compared to the non- 
linearity. However, a closer look at the adiabatic eigen- 
values and eigenstates shown in fig. [7] reveals a fatal non- 
linear avoided crossing scenario around t w —380, which 
will be referred to as the (avoided) "horn crossing" in the 
following because of its shape. In the linear case the sys- 
tem's dark state is rotated from ipds(t = — oo) = ipi to 
ipds{t = +co)) = "03, which leads to a coherent adiabatic 
population transfer. In the horn crossing however, this 
"dark state" disappears when it merges with a nonlinear 
eigenstate (which will be referred to as horn state in the 
following) , such that no adiabatic passage is possible any 
longer. To illustrate this breakdown of nonlinear STI- 
RAP, we integrate the three-level NLSE numerically for 
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FIG. 9: As fig. |3 however for g = +0.05. 

the same parameters as in fig. El The resulting evolution 
of the populations |a(t)| 2 and \c(t)\ 2 is shown in fig. in 
comparison to the population |a a d| 2 of the instantaneous 
eigenstates. Due to the crossing, a state initially pre- 
pared in level ipi, a(t = — oo) = 1, cannot be transfered 
to level ^3 adiabatically any more. The dynamics follows 
the instantaneous eigenstate (the "dark state") adiabat- 
ically until this state disappears at the horn crossing. 
Fast oscillations of the populations |a(£)| 2 and |c(£)| 2 are 
observed afterwards. Finally the system settles down to 
a steady state again, but the population has not been 
transferred completely. 

For g, A > 0, it is found that the nonlinear eigenstate 
(the horn state) that merges with the dark state is a 
superposition of ipi and tp2 alone in the limit t — > — oo 
and consequently v,w — > 0. Substitution of c = and 
v = w = into equation (@J immediately leads to 

-A + g\a\ 2 = f i = g\b\ 2 . (11) 

Solving for A yields a condition for the existence of the 
horn state, 

A = 5 (H 2 -|6| 2 )<5- (12) 

In this case a complete population transfer using the STI- 
RAP scheme will be prevented by the horn crossing see- 



6 




-600 -300 300 600 

t 



FIG. 10: As fig. IHJ however for g = -0.2. 

nario discussed above. Nonlinear STIRAP is still possible 
if the nonlincarity is smaller than the detuning, g < A. 
This is illustrated in fig. where the evolution of the 
populations |a(i)| 2 and |c(i)| 2 is plotted in comparison 
to the population |a ac i| 2 of the adiabatic eigenstates for 
g = 0.05 and A = 0.1. One observes that the dynam- 
ics closely follows the adiabatic eigenstate and that the 
population is transferred from ipi to ipz completely Note 
that the probability of Landau-Zener tunneling is also in- 
creased for g < A if the system parameters are varied at 
a finite velocity (cf. section HTT1 and ref. 0)- 

An analogous situation is found for <?, A < 0. However, 
a different crossing scenario arises if the signs of A and g 
are opposite. Again the "dark state" disappears when it 
merges with a nonlinear eigenstate (horn state). In this 
case it is found that the horn state is now a superposition 
of ipi and "03 alone in the limit t — > — 00 as the dark state 
and that it exists for all values of g. Rigorously speaking, 
this leads to a breakdown of STIRAP even for very small 
nonlinearities. However, this crossing turns out to be not 
as fatal as the one discussed above; the transfer is still 
close to unity for weak nonlinearities (cf. fig. . An 
example of the dynamics is shown in fig. 1101 for A = 0.1 
and g — —0.2. 

In conclusion we find the following conditions for the 
feasibility of a complete adiabatic population transfer us- 
ing the STIRAP scheme: 

3 A>0 and |<7|<<7 C =|A| (13) 

The dependence of the transfer efficiency \c(t — > +00) | 2 
on the nonlinearity is shown in fig. 1111 for A = 0.1 and 
the same couplings v(t) and w(t) used in fig. and [S] 
Note that exactly the same dependence is found if the 
signs of both A and g are altered. The transfer efficiency 
is slightly reduced for all values g < and shows an 



oscillatory behavior. For g > 0, one clearly observes an 
abrupt breakdown of the transfer efficiency above the 
critical nonlinearity g > g c = |A|. 
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FIG. 11: Numerically calculated transfer efficiency of nonlin- 
ear STIRAP in dependence of the nonlinearity g for A = 0.1. 
The arrows indicate the values of g for which the dynamics is 
shown in fig. 151 l9l and 1101 respectively. 



V. CONCLUSION 

In conclusion, the eigenstates and the dynamics 
generated by the nonlinear Hamiltonian J3J are analyzed 
for two important cases: the equal-slope Landau-Zener 
model and the STIRAP scheme. The emergence of new 
nonlinear eigenstates and novel crossing scenarios leads 
to a breakdown of adiabatic evolution if the nonlinearity 
\g\ exceeds a critical value. Consequently STIRAP fails 
if the nonlinearity exceeds a critical value given by the 
detuning or if the nonlinear parameter and detuning 
have different signs. A novel feature of nonlinear Zener 
tunneling compared to the two-level system is the oscil- 
latory behavior of the transition probability P(ce). Open 
problems include a detailed analysis of the oscillations 
of -P(a), the Landau-Zener scenario for non-equal slope 
and the effects of a parameter variation at finite velocity 
on the nonlinear STIRAP scheme. 
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